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Abstract 

Exact calculations are given for the Casimir energy for various fields in R x S" 3 geometry. The 
Green's function method naturally gives a result in a form convenient in the high-temperature 
limit, while the statistical mechanical approach gives a form convenient for low temperatures. The 
equivalence of these two representations is demonstrated. Some discrepancies with previous work 
are noted. In no case, even for M = 4 SUSY, is the ratio of entropy to energy found to be bounded. 
This deviation, however, occurs for low temperature, where the equilibrium approach may not be 
relevant. The same methods are used to calculate the energy and free energy for the TE modes in 
a half-Einstein universe bounded by a perfectly conducting 2-sphere. 
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I. INTRODUCTION 



The remarkable appearance of the holographic principle has fostered the understanding 
that some hitherto distant branches of theoretical physics may have a much deeper common 
origin than was expected. One outstanding example of this sort is the relation, suggested 
by Verlinde [l], between the Cardy entropy formula and the Friedmann equation for 
the evolution of the scale factor of the universe. Moreover, the proposal that there exists a 
holographic bound on the cosmological entropy associated with Casimir energy (for a general 
introduction to that subject, see Refs. |3|, f4|) suggests that there should be a deeper relation 
between Friedmann cosmology and the Casimir effect. 

There has much activity in the study of the Cardy- Verlinde formula and its various 
applications (see, for example, Refs. |J |6|, [7[ §, || and references therein). Specifically, there 
has been much interest in studying the entropy and energy arising from quantum and thermal 
fluctuations in conformal field theories 0, M. As we remarked, particularly interesting 
is the Verlinde bound |l] for the ratio of the entropy to the thermal energy. Whether this 
bound can be realized in realistic situations is a matter for specific calculations. Previous 
computations [§, |9| have been limited to the specific regime of high temperatures, so they 
are unable to provide definitive results. Here we obtain exact results for various fields in 
the Rx S 3 geometry, so the issues may be more decisively addressed. We also will consider 
electromagnetic modes in a half-Einstein universe, which give similar results using the same 
methods. 

The outline of this paper is as follows. In the following section we consider a conformally 
coupled scalar in the Rx S 3 geometry. We obtain exact results for the Casimir energy and 
free energy using both Green's function and statistical-mechanical methods. The former 
result is most convenient for the high-temperature regime, while the latter is suited to low 
temperatures. However, the results obtained by the two methods are exactly equal, as may 
be shown either by the use the the Euler-Maclaurin sum formula, or by the Poisson sum 
formula. For the physical applications to the evolution of the universe, except, perhaps, for 
the very earliest stages, the high-temperature results are the most relevant. In Sec. HI], for 
comparison, we compute the same results for R x S 1 . Vectors and fermions are treated in 
Sees. [TV] and [V], respectively. Entropy bounds are the subject of Sec. In the presum- 
ably physically relevant high-temperature regime, the ratio of entropy to energy is small, 
of order (aT)" 1 , where T is the temperature and a is the radius of the sphere. Linear 
temperature terms in the energy, omitted in earlier analyses [||, preclude the existence of 
an entropy bound even in the high-temperature regime. At low temperature, moreover, 
both the entropy and the thermal deviation of the energy from the zero-point energy be- 
come exponentially small, while the ratio of these two quantities diverges as T — > 0. In 
Sec. [VII| we consider the half-Einstein universe. Some summarizing remarks are offered in 
the Conclusions. There are three Appendices. The first points out the fundamental diffi- 
culty of breaking conformal symmetry in such problems, the second shows the applicability 
of the Euler-Maclaurin method to the example treated in Sec. [VI I , and the third responds 
to a suggestion that in R x S 1 the zero-mode contributes a linear temperature term at low 
temperature. 
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II. CONFORMALLY COUPLED SCALAR 



A. Green's Function Method 



We can start from the formalism given in Kantowski and Milton |TIJ, which in turn was 
based on the Green's function formulation first given by Schwinger et al. [IT|. The energy 
is given by the imaginary part of the Green's function, 



U = V 3 d d'°SG(x,y;x',y' 



i x=x' ,y=y 



(2.1) 



where the "external" coordinates x consist only of the time. We introduce a Fourier trans- 
form there 

lLC ' ' ' ), (2.2) 



G(t,y;t',y' 



271 



1 ~'g(y,y';v), 



so 



(2.3) 



where the contour c encircles the poles on the positive axis in a negative sense, and those 
on the negative axis in a positive sense. The reduced Green's function satisfies 



(Vg + uj 2 )g(y, y'; u) = -5(y - y'). 
The eigenvectors and eigenvalues of the Laplacian on S 3 are 

M? 



V^ m (|/) 



-Y t ™, Mf = 1(1 + 2), 



which have degeneracy D\ = (I + l) 2 . The addition theorem is 

J2 Y r(y) Y r(y) 

m 

so in view of the eigenfunction construction 



El 

v 3 



(2.4) 



(2.5) 



(2.6) 



, ^ Yr(y)Yr(y') 



lm 



the Casimir energy is 



U = — 



Air 



duu 2 



D, 



Mf/a? - u 2 



(2.7) 



(2.8) 



Conformal coupling in incorporated by replacing the Laplacian in Eq. (2^) by V| + £R, 
which for the conformal value £ = 1/6 amounts to the addition of 1 to the Mf operator, 
that is, replacing Mf —>•(? + I) 2 - 

Temperature dependence is incorporated by the replacement 



du 4i 
2^ ~* /? 



E'. 

n=0 



2rm 



(2.9) 



The prime on the summation sign means that the n = term is counted with half weight. 
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We carry out the sum on / in Eq. fl2.8|) as follows: The general representation 

oo 

El 7T 1 . 
— - = cotTra (2.10) 

becomes here 



m 2 — a 2 2a 2a 1 

m=0 



E(l + l) 2 7T . „ 7lUJa 3 t . 
- — \ n , ' = constant H — iioa -\ — . (2.11) 
+ 1 ) 2 a 2 -u 2 2 e -2^a _i y ) 

1=0 K 11 

We then make the finite-temperature replacements ( p.9[ ) and obtain, after dropping the 
contact term arising from the constant (oc £(— 2) = 0) 



U 



1 /27ra x ; 



a V P 



V - + — 

e Air 2 an/P _ J 240 

n=0 



(2.12) 



which, since the summand vanishes at n — 0, gives only exponentially small corrections to 
Stefan's law, 

rr 1 (2vraT) 4 



Here, and throughout the paper, we use the formulas [T^ 



C(l-2n) = - — , C(2n) = |5 2n L n = 1,2,3,..., (2.14) 

sv ) 2n sv ; 2(2n)! 1 1 ' v ; 

i?2n being the Bernoulli numbers. 

B. Statistical-Mechanical Approach 

We recall the usual statistical mechanical expression for the free energy, 

F = -kTlnZ, (2.15a) 

InZ = - ln (! - e_/?Pl ) 
j 

oo 

= _ + 1)^2 ln (1 _ e -/J(n+l)/a) ; ( 2 15b ) 

n=0 

for conformally coupled scalars in Here the zero-point energy has been subtracted. 

The latter, of course, is easily calculated, 

* = + !)<->=+_ = -C(i- rf ). (2.16) 

n=0 

The specific results for two and four dimensions are 



For the temperature dependence, we differentiate the partition function, 



ly 2 f , 5= A. (2.18) 

n=l 



d/3 a 

This is a very different representation from Eq. ( |2.12| ). Nevertheless, from it we may obtain 
the same result we found above if we use the Euler-Maclaurin sum formula, 



x POO 1 

£)/(n) = / dnf(n) + -[f(oo) + f(0)} 

n=0 ^° 



+ £ [/^W - / (2fc - 1} (0)] • (2.19) 



Doing so here yields for the integral 



a jo 

while the third derivative term gives 



1 m , n 3 a 3 vr 4 

dn-s-. = — ~~r , 2.20 



3B 



4 



a 4! 240a ' 



(2.21) 



so we exactly reproduce the results of Kutasov and Larsen || , but derived by a much more 
transparent method, 

' 1 



coinciding with Eqs. Q2.13Q and (|2.17b| ). 



C. Comparison with Usual Finite-Temperature Casimir Effect 

The same argument could be used for a scalar field subject to Dirichlet boundary condi- 
tions on parallel plates, with d transverse dimensions (see Ref. j3J): 

F = TV f J?*I ± In (l - e -V^») . (2 . 23) 

^ ' n=—oo 

If we use the E-M formula, and replace the sum by an integral, and rescale variables, j3k — > k, 
f3im/a — > n, we obtain 

F= — [*t^Ul-e-^) 

~ P d ^ J {27TY+ 1 V J 
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where A d = 2n d / 2 /T(d/2) is the area of a d-dimensional sphere. This result is not the correct 
linear high-temperature T dependence, but instead its negative appears as a (subdominant) 
term in low-temperature limit! (See, for example, Ref. Eq. (2.91).) What is different 
about this case is that there is an inside (between the plates) and an outside, and the Stefan's 
law term, being proportional to the volume, cancels between the two regions — it can give 
no force on the plates. Here there is no outside, so the volume term must be included. As 
explained in Ref. ||], p. 56, the high temperature dependence for parallel plates of area A is 
linear, 

F _ Cg) 
A ~~ 167ra 2 
while the low-temperature correction is cubic, 



d = 2 : — = -T-^-^, aT > 1, (2.25a) 



d = 2: ^ = -T^-^ + ^T\ aT«l. (2.25b) 



A 1440a 3 4vr 90 

(Note that the cubic term does not contribute to the force on the plates, being independent 

KM 

4-7T 



of a, but does contribute to the entropy, S ~ Kl 3 !^ 2 ^ a T <^ 1.) 



D. Relation Between Representations 

We have two representation for the Casimir energy, the one obtained from the Green's 
function, Eq. ( 2.12 ), and the one obtained from the partition function, Eq. (|2 .18 ). The rela- 
tion between the two can be found from the Poisson sum formula. If the Fourier transform 
of a function b(x) is defined by 



cl a 



/ —e-^bix), (2.26) 

J-oo 27r 



then 

oo oo 

K n ) = 2?r Yl c ( 2nn )- ( 2 - 27 ) 

n=— oo n=— oo 

This is just a rewriting of the elementary identity 

oo oo 

£ e - i2nnx = S(x-n). (2.28) 

n=— oo n=— oo 

So start from Eq. ( p. 18 ) and consider 

1 f of T > n 

b( x ) = - \ e'^^-v x ^ u ' (2.29) 



a { 0, x < 0. 
It is easily seen that the Fourier transform is 

i L I i ) L / a \ ' 

c(a) 



i r 4 i r a y /JX>/ . . 

— ? tttti rr tj = t; r(4)C (4, 1 + iaaf/3) , 2.30 

2vra^ \p (k + 1 ) a + ia] A 2ira \fi J VAV ' /M; ' 
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in terms of the Hurwitz zeta function. Thus the energy is 



e = - ( % ) r(4) y y - 4 . (2.3i) 



If we sum this on n first, we obtain the alternative expression 



E=h{A)-±- v A y C(4,l-W(2vra)) + C(4) 

° \ 71 ' lfc=-oo 



27ra x 4 



/3 



(2.32) 



By comparing with the original expression ([2.18 ) for E we can replace the sum over Hurwitz 
zeta functions by the expression ( [2.12] ): 

1/2™ \ y^L. — + J_ r (27raT)4 _ x i (2 33) 

x ' n=l 

to which must be added the zero-point energy, which cancels the T-independent term. (Thus 
the T = term, which is the ZPE, arises from the sum, again as directly verified by the 
E-M formula, and in the low temperature limit the Stefan's law term is cancelled.) Thus 
we again reproduce Eq. Q2.12 ). The two representations are adapted explicitly for the two 



limits: 

U=— + ±y—^ (2.34a) 

240a a ^ e 2 ™ 5 - 1 v ; 

71=1 

1 1 00 3 

= (27raT) 4 + -(27raT) 4 y — \ . (2.34b) 

240a v ' a y ' e 2 ™ a / 5 - 1 v ' 

n=l 

III. d = 2 CONFORMAL SCALAR 

Because there is a subtle issue involving zero-modes here, it is useful to repeat the above 
calculation for d = 2. The partition function is given by 

00 

lnZ = -^ln(l-e^ (n+1)/a ). (3.1) 

n=0 

This immediately gives the low-temperature representation, 

I 1 00 

II n 

U = -^ + a^e 2 ™ s -V (3 ' 2) 

71=1 

displaying an exponentially small correction to the zero-point energy if (3 ^> 1. Again, by 
use of the Euler-Maclaurin sum formula we can obtain the high-temperature limit, 

U ~ i (27raT)2 "^ T ' aT>>L (3 ' 3) 

This again coincides with the result found in Ref. ||. However, the linear term in T is 
omitted in the analysis of Klemm et al. |9j because they believe that the result can only 
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be trusted in the limit T — > oo |13 |. This omission is necessary for their derivation of the 
Cardy formula [Q]. 

It is easy to reproduce this result from the Green's function method. After the finite- 
temperature substitutions, the expression is 



% 8iir ^—^ , f 2im\ air 



' n=0 v r ' 



1 + 



g47r 2 an//3 _ ^ 

1 1 1 n 

24^P~2 "^^e 2 ™/*-!' ( ' J) ' ' 



n=l 



which gives the explicit exponential corrections to the high temperature limit ( |3.3j ). The 
low-temperature limit displayed in Eq. ( p.2| ) may be easily obtained from this by using the 
Euler-Maclaurin sum formula. 

A proof of the equivalence of the two representations (|3.4j) and (|3.2j ) can be carried out 



along the lines sketched in Sec. [II D| . There is one complication, due to the presence of the 
zero-mode at n = 0. The Poisson sum formula (|2.27|) , applied to a function b(x) which is 
positive for x > and zero for x < implies a Fourier transform c(a) which has poles in 
the upper half plane and hence implies b(x) = for x < 0. This contradiction at x = is 
resolved by the following prescription: Write 

OO OO 

5>(n)=X>(n)--6(0), (3 ' 5) 

n=l n=0 

replace the latter sum by the sum over the Fourier transform c(2nn), and then when we 
interchange the sums and transform back the zero mode is not present. Thus only the term 
explicitly displayed above arises from the n = term. In this way the equivalence between 
the two representations (|3.4|) and ( |3.2|) , with the —\T term, is established. 

IV. VECTOR FIELD 

The analysis proceeds similarly to that given in Sec. II. For S^ -1 the degeneracy and 
eigenvalues are 

21 ft+f-1) (l + d-2)(l + d-4)\ 
Dl ~ (d-3)(Z + l)! ' 

Mf = 1(1 + d- 2), (4.2) 

so for d = 4 if we again add the conformal coupling value 1 to Mf we obtain for the Green's 
function mode sum 

21(1 + 2) 1 2 / 1 \ / 2 \ 

> 71 xo / n n ~> o + ma H n~ • ( 4 -3) 

(I + I) 2 /a? -uj 2 u 2 V uoa) \ e ~ 2mwa -lj K 1 

After making the finite temperature replacement, we carry out the sum on n, with the result 

(2naTf 1 (2ixaT) 2 



U 



120a 12 a 

2(2naT) 2 ^ r , v9 1 
+ T + ~^~ n E [" + (2-aT) V] pAM _ r (4.4) 



n=l 
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where the T term comes from the n = term in the sum. [Note that the leading term is 
exactly twice that found for a scalar field found in Eq. ( j2.13|) , as we would expect.] Since 
the remaining sum is exponentially small in the large T limit, this form is well-adapted for 
high temperature. (The T 4 and T 2 terms are as given in Ref. ||, while the linear term is 
implicitly given there.) However, it is exact, and by using the Euler-Maclaurin sum formula 
it yields the low temperature limit, 

u ~T^< aT<<1 ' < 4 ' B) 

up to exponentially small corrections. The latter may be directly inferred from the partition 
function, 

CO 

In Z = - 2l ( l + 2 ) ln i 1 ~ e" w+1)/a ) . (4.6) 
i=i 

By taking the negative derivative of this with respect to (3 we obtain the alternative repre- 
sentation 

u = JL + l±%zJl. (4 . 7) 

120a a ^ e^/ a - 1 V 1 

i=i 

The Euler-Maclaurin formula applied to this last sum yields the leading high-temperature 
term in Eq. (fOJ), including the term linear in T, and the exact equivalence of the two 
expressions ( |4.7[ ) and ( |4.4| ) again is demonstrated by the Poisson sum formula, including the 
zero-mode prescription enunciated in Sec. |TTT|. 



V. WEYL FERMIONS 



The procedure is now routine, except for the complications due to fermions. The degen- 
eracies and eigenvalues are 



A = 2(1 + 2)(l + 1), M? = 1(1 + 3), 



(5.1) 



so recalling the minus sign associated with a fermionic trace, and that the fermionic thermal 
Green's functions must be antiperiodic, we have the following expression for the energy, 



U 



l o7TZ 



4vr (3 



\- ^ ma 

2a 2 > 

^ 2 



n=0 
X ( 1 + 



i2ix(n + 1/2) 



1 /z2tt(™ + 1/2) 
3 V P 



8vr 2 a(n+l/2)//3 



= 47r 2 a(n+l/2)//3 



.r 4 - -r 2 

a 960 96 



1 °° / 9 

— E K 2 - + ^ + (2- + l)^ 1 ] eM2n+ iys _ i - 

n=0 ^ 



3 vr(2n+l)/5 



• (5-2) 



The low-temperature limit (the zero-point energy) may be obtained from this by the Euler- 
Maclaurin formula, and the exponential corrections in that limit may be obtained directly 
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from the partition function, 



In Z = 2 < n + !) ln i 1 + e- /3(2n+1)/2a ) . (5.3) 

That is 



n=l 



U 



1 

a 



17 ^n(n + l)(2n + l) 



(5.4) 



960 ^ g/3(2n+l)/2a _|_ J 
n=l 

The equivalence between Eqs. ( |5.2|) and (|5.4|) may be demonstrated again either though the 
Euler-Maclaurin sum formula applied to the latter, or exactly through the use of the Poisson 
sum formula. 



VI. ENTROPY BOUNDS 

From the above results, thermodynamic information may extracted in terms of the free 
energy (with zero-point energy E subtracted), 

F = -T\nZ, (6.1) 

in terms of which the energy 

E S U-K—±h>Z=pF, (6.2) 

and the entropy 

S = 2na5 2 ^-F = f3(E - F) (6.3) 
do 

may be extracted. 

A. Two-dimensional scalar 

Klemm et al. || ignore the linear T term in the energy, and so have for 5^1 (see 



Sec. imp 



E = -L(«T 2 + 1), (6.4a) 



24a 



24a 



1 (r 2 -i), (6.4b) 



6 

These imply the Verlinde-Cardy formula [[[], [| 



S = -5-\ (6.4c) 



S = 4ira^/-E (E + E ), (6.5) 

and the entropy bound 

2-^- < 1. (6.6) 



2-naE 5 2 + 1 
10 



However, this result is not meaningful as it stands. Even in the high-temperature limit we 
must add the term linear in temperature to the energy, which implies instead from Eq. Q3.3| ) , 
for 8 < 1, that 



S = ^r' + ^lni-l). (6.7c) 

The ratio of S 1 to E is then unbounded as 5 — > oo. Yet this takes us to the low-temperature 
regime, where we must use the leading exponential corrections, for 8 ^> 1, 



£ ~ -e"^, (6.8a) 
a 

F 4 e-/3/a > ( 6 - 8b ) 

5 „ t e -M* % (6.8c) 



so the entropy-energy ratio is 



S 



2,aE = «• «> L < 6 ' 9 » 
It is apparent that this latter result is universal because the energy always dominates the 
free energy in the low temperature regime. 



B. Entropy Bounds in Four Dimensions 

In the following we will consider cases with N s conformal scalars, vectors, and Nf 
Weyl fermions. For example, the M = 4 supersymmetric multiplet has N s = 6, N v = 1, and 
JV) = 4. 

In the high-temperature regime we may write the free energy, energy, and entropy as 1 

F ~ --\a A 5~ A + a 2 8~ 2 + a^lnS + a ], (6.10a) 
a 

E ~ ^[3a 4 r 4 + a 2 r 2 - a - aid' 1 }, (6.10b) 

5 ~ 27r[4a 4 r 3 + 2a 2 r x - a x {l - \n8)}. (6.10c) 

Here the coefficients were determined in the previous sections to be [see Eqs. ( |2.33|) , ( |4.4D , 



What is called the Cardy formula is simply the observation that the leading behavior of S is the geometric 
mean of the leading and subleading terms in E. The term "Casimir energy" for the latter is misleading 
in other than 1 + 1 dimensions. The entire energy U is due to quantum and thermal fluctuations, so it 
all should properly be reckoned as Casimir energy. 
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and Q] 

N N 7N f 

cla = — , 6.11a) 

720 360 2880' v ; 

a2 = -12-%' (6 ' llb) 

ao = 3a4 — 02, (6.11c) 

iV 

ax = (6.11d) 

Even ignoring the a\ term, Klemm et al. || note that no entropy bound is possible, unless 
special choices are made for the field multiplicities. For the M = 4 case the first three 
coefficients are 

1 13 

dd = — , a,2 — — j a o — — i (6.12) 

48 8 16 v ; 

and the entropy-energy ratio becomes 



S 1 -ln5 + fr 3 (l-35 2 ) 

2^E ~ 5- 1 + |5- 4 (l + 5 2 )(l-35 2 ) 



(6.13) 



If the ai terms here were omitted, the zero in both the energy and entropy at 5 2 = 1/3 
would cancel, and we would have the limit given in Ref. 0: 

5 4 4 <i (6.14) 



2naE 3 1 + 5 2 ~ 3 

in the high temperature regime. But a\ ^ 0, and the ratio ( 6.13|) diverges as 5 — > oo. Of 



course that limit is the low-temperature one, but the argument given in Sec. |VI A| then 
applies and shows that 

' S 5, 5^oo. (6.15) 



2naE 

Although in this limit both the entropy and the subtracted energy are exponentially small, 
their ratio is unbounded. 

The reader should note that we are not in formal disagreement with previous studies 
|8|, [|. The interest there was restricted to high temperature, which is presumably all that 
is relevant to nearly the entire history of the universe. 2 In that case, only the leading terms 
in 1/8 are relevant, and the ratio of entropy to energy is always of order i<l. It is not 
surprising that such results as Eq. (|6.14D are an unreliable guide to the moderate and low 
temperature regimes, which might be relevant in the very earliest (pre-inflationary) stages 
of the universe. 

Another point, which is more closely connected with physics, is that it is permissible 
to make use of the thermodynamical formalism for fluctuating quasi-classical systems only 



2 It might be noticed that before or after photon decoupling, but after inflation, the value of aT stays nearly 
constant Jlij , essentially reflecting entropy conservation. That value is far into the high temperature 
regime, the present value of S being Sq ~ 10~ 30 . Insofar as it is permissible to speak of temperature 
during inflation, aT is also constant then, but of a much smaller value, which value increases dramatically 
during reheating. 
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when the temperature T is sufficiently high. As discussed recently by Das et al. ||15|| , their 
Eq. (12), one must have 

T>~, (6.16) 
r 

where r is the relaxation time, in order to use thermodynamics. Now, in our case it is 
presumably legitimate to estimate r to be of the same order of magnitude as the transit 
time for light across a distance of order a, i.e., r ~ a. This leads to the condition T> 1/a, 
which actually means 

<J<1. (6.17) 

In other words, the physical condition for using Eq. (|6.10 ) above seems to be that the linear 
term, which is of order 5 3 relative to the first term, is negligible! It seems that we must be 
careful in not assigning too much physical significance to the subleading corrections. 



VII. ON THE HALF EINSTEIN UNIVERSE 

The metric of the static Einstein universe has continuously attracted interest, both be- 
cause it is easily tractable analytically and also because it is conformally equivalent to all 
the closed Robert son- Walker metrics. The Einstein metric can be written as 

ds 2 = -dt 2 + a 2 [d X 2 + sin 2 X (d0 2 + sin 2 9 d(j)% (7.1) 

where 9 G [0,7r] and G [0, 27r]. In the case of the full Einstein universe, x [0,vr]. The 
energy density p consists of two parts, one a matter (dust) part, p , and one a vacuum 
part, A/8nG, A being the cosmological constant. The pressure is p = —A/8ttG. From the 
Friedmann equations, po = A/47rG, and the scale factor becomes a = hr l l 2 = (AirGpo) ~ 1 ^ 2 . 

The half Einstein universe, in which x varies only from to tt/2, turns out to be an 
interesting variant of the Einstein static universe idea; see, for instance, Refs. and |Ty] . 



As discussed in Ref. ||18|| , we can consider this universe as a three-dimensional spherical 
volume, spanned by the "radius" x an d the angular coordinates 9 and 4>, closed by a two- 
dimensional spherical surface lying at x — 7r /2- Let us assume that this surface is perfectly 
conducting. This is the simplest imaginable option in the electromagnetic case, and is 
analogous to the Dirichlet boundary condition in scalar field theory. The assumption allows 
the consideration of standing electromagnetic waves in this spherical cavity. We expect to 
find the same possibility of dividing the possible eigenmodes into independent TE modes and 
TM modes, as we do when dealing with the conventional standing modes in the Minkowski 
metric. Explicit calculation of the fundamental electromagnetic field modes was actually 
given in Ref. Jl8| . 

Let us assume an orthonormal basis, designated by carets, {u 1 , u*, u 9 , u^} = 
{dt,adx, asinxd9, a sin x sin 9 d(p}, and split off the time factor as e~ luJt . From Maxwell's 
equations in the curvilinear space we find the governing equation for E% (or H%). Denot- 
ing these field components collectively by X, and writing E%(x, 9, (ft) = i?^(x)Y' z m (^, (f>), we 
obtain 

d 2 

— (sin 2 x X) + (ujo) 2 sin 2 X X - 1(1 + 1)X = 0. (7.2) 
dx 

The solution of this equation is known: 

Xocsin^xC^^cosx), (7.3) 
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where n is an integer, being the Gegenbauer polynomials [12|. The differential 

equation satisfied by Cp a \x),p > 1 is 

(1 - x 2 ) C^'\x) - {2a + l)xC^\x) + p(p + 2a) C { p a \x) = 0. (7.4) 
Inserting Eq. ( |7.3| ) into Eq. ( |7.2| ) we obtain the eigenfrequencies 

n + 1 

u n = , n > I > 1. (7.5) 

a 

For the electromagnetic field, we know that I > 1. Moreover, in order to avoid infinities at 
the origin x = 0, we must have n — I > in Eq. (|7.3|) [[18] . 



As mentioned above, the TE and TM modes are independent and can be considered 
separately. In the following, we consider the TE modes only. 



A. The TE Modes 

The electromagnetic boundary condition at the conducting surface is in this case 

Hx = 0, X = tt/2. (7-6) 

Since / > 1 the condition ( [T.6| ), together with the general property that Cp (0) = when p 
is odd fl8|| , shows that the subscript {n — I) in Eq. ( |7.3j ) must be odd. The eigenvalues ( |7.5|) 



are seen to depend only on n, so that they are degenerate with respect to / as well as to the 
magnetic quantum number m appearing in Y^ m . We accordingly have to sum (21 + 1) over 
all the admissible values of /. 

Assume first that n is even. According to the above, / can then take all odd values 
between 1 and (n — 1), so that the degeneracy becomes 

(2Z + l) = ^(n+l), 71 = 2,4,6,.... (7.7a) 

;=l,3,5,..(n-l) 

Next, if n is odd, I can take all even values between 2 and (n — 1). This leads to the 
degeneracy 

]T (2/ + l) = !^zl(n + 2) ) n = l,3,5,.... (7.7b) 

Z=2,4,6,..(n-1) 

The partition function for the TE modes is accordingly 

oo 

In Z TE = - ~( n + !) ln i 1 ~ e^ (n+1)/a ) 



2 

n=2,4,6,.. 



J2 ^(ri + 2) In (l-e-^ +1 )/ a ; 



n — 1 

n=l,3,5, 
oo 

- n ( 2n + !) ln t 1 - e" /3(2n+1)/a ) 



n=l 

oo 



^(n-l)(2n + l)ln(l-e- 2/3n/a ). (7.8) 



n=l 
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We write this expression as a sum of two terms, 

lnZ TE = A + B, (7.9) 

where A, the first term in Eq. ( |7.8| ) , is rewritten as 

-.00 , 00 

A = -- ^(2n + l) 2 In (l - e -M*»+U*) + i J^(2n + 1) In (l - e - 27r ( 2n+1 ) 5 ) , (7.10) 



/ 2 

n=l n=l 



where we recall 5 = f3/{2na). The generic form is thus 

00 

= ^(2n + l) a ln(l-e^ 2n+1 ) <5 ), a = 2,1. (7.11) 



n=l 

For variety's sake, we evaluate this expression using a variation on the method given 
in Ref. ||. (We also used this method to derive the high temperature results in Sec. II.) 
Expanding the logarithm and thereafter differentiating with respect to 5, we have 

= 2vr^(2n+ l)" +1 ^e- 2 ^ 2 ™ +1 ) 5 , (7.12) 

n=l k=l 

which can be processed further by inserting 

e~ x = — ! dsx- s T{s). (7.13) 

The integration is here taken along a line parallel to the imaginary axis with sufficiently large 
value of dts. By introducing Riemann's zeta function ((s) via the sum over k we obtain 



dG^ _ 1 
85 i 

The sum over n leads to 



00 „ 

V(2n+1) Q+1 / ds[2Tt(2n + l)5}~ s T(s)((s). (7.14) 

n=l ^ 



^(2n + l)-*-*"* 1 = (1 - 2- s+Q+1 )C(s - a - 1) - 1, (7.15) 



n=l 



and the expression ( |7.14| ) can conveniently be written as 



05 
with 



' = {I}<«> + {!!}, (7.16) 



{/}(«) = I / ds(27i5)- s (l-2- s+a+1 )((s-a-l)T(s)((s), (7.17a) 
z ic 

{//} = -~ I ds{2it8)- s T(s)({s). (7.17b) 
1 Jc 
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Consider first the expression ( |7.17a|) . Because of the poles of the zeta functions it is properly 
defined only for !fts>a + 2(a>0 assumed). For lower values of s, we define the zeta 
functions in terms of their analytic continuations. The integrand has poles at s = 1 and 
s = a + 2 from the zeta functions, and at s = —n, n = 0, 1, 2, ... from the gamma function. 
There are, however, zeroes from the zeta functions, at s = —2n and s = a + 1 — 2n, 
n — 1,2, 3, so that the remaining number of poles becomes effectively reduced. 

For a = 2, there remain poles only at s = 4 and s = 0. We distort the contour to encircle 
these poles and find, by means of the formulas fl2.14| ), that 

n\m = -^ + —. (7.18) 
1 1 2405 4 120 K ' 

For a — 1, the expression ( pM7a| ) has poles at s — 3, 1, —1, —3, .... Only s = 3, 1 need to be 
taken into account since we assume that S -C 1 (high temperatures), implying that terms of 
order 5 are negligible. We get 

Consider now the expression {//}, Eq. ( 7.17b|) , which is independent of a. The integrand 
has poles at s = 1,-1, —3, —5, .... Because of the smallness of 5 we have to include only the 
lowest value, s — 1. We get 

{//} = -i (7.20) 

It actually turns out that when constructing OA/ '85 according to Eqs. ( |7.10| ), ( |7.11[ ), and 
( [T.16| ) the quantity {//} does not contribute: 

d_A = _l_dGM IdGM = _1 (2) 1 (1) 

35 2 85 2 05 2 X s 2 X s 

t C(3) 1 7tt /r70iN 
+ + 7771 - (7-21) 



4805 4 8tt 2 5 3 245 240 
We now consider the second term B in Eq. (|7.9|): 



B = - J^(2n 2 - n - 1) In (1 - e - 47rn<5 ) , (7.22) 

n=l 

whose generic form is 

oo 

#(°) = ^n Q ln(l-e- w ) , a = 2, 1,0. (7.23) 

n=l 

A calculation along the same lines as above yields 



- / ds{A7i5)- s C{s-a-l)T{s)C{s), (7.24) 



05 

and we get (the values of s giving poles are written within the parentheses) 



OH^ 


71 


71 


05 


960 5 4 


~ 60 


8HW 


C(3) 


1 


05 


8vr 2 5 3 


12 5 


dH (0) 


71 


1 


05 


24 5 2 


-2S + 



(s = 4,0), (7.25a) 
(s = 3,l), (7.25b) 
I (s = 2,1,0). (7.25c) 
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Thus 

OS 480 5 4 8tt 2 5 3 24 5 2 125 5' y ' ' 

Adding Eqs. ( |7.21| ) and ( |7.2(j| ) we thus get for the thermodynamic energy E TE of the TE 
modes 

„ TP 3 , 71 C(3) 7T 13 417T ,„ „„ s 

2 ™ £ = -« 2 = 240F - - 2«5 + - 240 ■ (7 ' 27) 

This corresponds to the free energy 

The same result is obtained in Appendix [B] using the Euler-Maclaurin method explained in 
Sec. II. 



B. Comment on TM Modes 

Although we have not considered the TM modes, we expect on physical grounds that 
they should contribute the same amount to the leading term in the energy as do the TE 
modes. Accordingly, for the total thermodynamic energy U = U^+U™ in the half Einstein 
universe we expect 

7T 

27iaU = — — — — rr + subleading terms. (7.29) 
120 o 4 

This is the same leading term expression as obtained for a conformal scalar field in the 
fix S 3 geometry considered in Sec. II, and half that for a vector, see Eq. ( |4.4|) . We conclude 
that in the high-temperature limit, the ratio of entropy to energy is the same as given in 
Eq. (EH , 

SSb-I* S<<1 - (7 - 30) 



VIII. CONCLUSION 

There are two main points in the Verlinde proposal PJ. First, the Friedmann- Robertson- 
Walker (FRW) equation may be written formally very similarly to the Cardy entropy equa- 
tion 0. Even in the presence of higher derivative gravitational terms, or unusual matter 
content, or in the presence of quantum corrections, one can redefine the energy density to 
arrive at the standard form of the FRW equations. 

Second, the free energy for higher dimensional conformal field theory has formally (at least 
for high temperature) the same structure as in the two-dimensional Cardy case. Moreover, 
assuming that (a closed) universe has extensive and subextensive ( "Casimir" ) contributions 
to its energy and entropy it was suspected that the "Casimir energy" is less than the 
Bekenstein-Hawking energy. As a result, the entropy/energy cosmological bound appears. 
It is, however, clear that this bound cannot be universal. For example, quantum corrections 
or the choice of an unusual state of matter would significantly alter such a bound. 

Our main qualitative result in this paper is that entropy/energy bounds should be relevant 
only in the ultra-high temperature limit, which applies to the universe after inflation. With 
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the decrease of temperature the bound becomes much less reliable, until at low temperature 
(aT <C 1), which might be the case in the very early universe, the entropy dominates the 
energy. This effect occurs already for conformal matter. As is apparent from the results 
of Appendix |A] the situation for non-conformal matter is much more complicated. Hence, 
it is unclear if entropy/energy bounds should exist at all even for high temperature. This 
question will be discussed elsewhere. 
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APPENDIX A: CONFORMAL SYMMETRY BREAKING 

In this Appendix, we wish to consider conformal symmetry breaking by choosing for the 
eigenvalue of the scalar field 

M / 2 = /(/ + 2) + X =(/ + l) 2 + X -l. (Al) 

The conformal value considered Sec. II corresponded to \ = 1. As we shall see, the analysis 
is considerably more subtle when \ ^ 1. 

Let us first attempt to calculate the zero-point energy, the Casimir energy at zero tem- 
perature. Using Eq. (|2.8|) and carrying out the sum according to Eq. (|2.10|) , we have, after 
dropping the constant, 

Eq — — j duu 2 a 2 a cot na, (A2) 



where 

a = ^u 2 a 2 + 1 - x- (A3) 

Remember that c encloses the poles on the positive real axis in a clockwise sense, and those 
on the negative real axis in a counterclockwise sense. There is no pole at a = 0, but poles 
at a = ±1, ±2, .... Changing to the variable a, we have 

E " = i I aHa Vq ' 2 + * - 1 (p^in + 1) ■ < A4 > 

Here we draw branch lines on the positive and negative real axes, starting from the branch 
points at a = ±-y/l — x, respectively, and d encircles those branch lines in a =F sense 
respectively. 

Now we see a serious problem. When % = 1 there is no branch line, and the 1/2 term 
goes away, as it is not singular. But with x ^ 1 each individual term possesses a branch line 
(although the entire integrand in Eq. ( |A4| ) does not), and this expression seems irreducibly 
divergent. 

Nevertheless, let us press on a bit. The finite temperature version of this is 

00 /r> \ 2 



= £ v ' ( 2nan \ 
' P 2 -! v J 



a cot Txct 




27r A /(27ma//3) 2 _ 



(A5) 



18 



Now we have further difficulties. What is one to make of the leading divergent term? Also, 
the square root is not real if a/ (3 <C 1. Ambiguities also abound if we start from the 
partition function, and try to derive from it the high-temperature limit. It seems apparent 
that conformal symmetry breaking poses problems for Casimir calculations that have not 
yet been solved. This is not a new observation, of course. For example, it is known [2C] that 
a massive field in flat space bounded by a spherical shell possesses divergences that are not 
present for massless fields. 



APPENDIX B: EULER-MACLAURIN METHOD 



In this Appendix we sketch the derivation of the result (7.27) using the Euler-Maclaurin 
method given in Sec. II. We write the result for the subtracted energy 



E 



a i 00 

-|^ TE = -$>W+SW], 



n=l 



where 



fin) 



n(2n+ If 

e /3{2n+l)/a _ I ' 

n(2n + 3)2(n + 1) 

e 2(n+l)/3/a _ \ 



(Bl) 

(B2) 
(B3) 



To get the high-temperature limit, we use the Euler-Maclaurin sum formula. The first 
integral term is 



e u - 1 



(B4) 



For /3/a <S 1 this is easily seen to be 

7T 

dnf(n) = — ( - 



4 / \ 4 

7r / a 



f3/a 

C(3) fa^ 3 

J, 



1 a 1 
Mp ~ 96' 



60 \(3J 2 

In the same way, the integral of the second function is seen to be 

, . , vr 4 faV C(3) fa\ 3 7T 2 fa\ 2 5a 1 

The contribution of these two integrals to the energy is 

i 



(B5) 



(B6) 



ir lay f a\ 3 vr 2 / a\ 2 7 f a\ 11 

30 W ~ W ~12W + 8\PJ~ 32 



(B7) 



The derivative terms in the Euler-Maclaurin formula may be obtained from the small j3 
forms: 



f( n ) = ^ n (2n + 1) - \n(2n + l) 2 + O (S- 
p 2 \a 

g{n) = ^ n (2n + 3) - n(2n + 3)(n + 1) + O ^ ) . 



(B8) 
(B9) 
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Thus, for high temperature, the only nonzero derivatives appearing are 



/(0) = 0, /'(0) 



p 



1 

2" 



f"'(0) 



-12, 



g(0) = 0, g'(0) = 3^-3, </"(0) = -12. 



(BIO) 
(Bll) 



(-24) 



These give a contribution to the energy of 

b a [ 2 V P V 24 
If a 31 
" a \~3j3 + 120 

Combining this with E a , Eq. ([B7|) , we obtain exactly the result found in Sec. VII, Eq. ( [7.271) : 

1 C(3) 1 13 41 



(B12) 



2S=i 
a 



4805 4 8vr 3 5 3 485 2 48vr5 480 



(B13) 



APPENDIX C: ZERO-MODE CONTRIBUTION 

There are two approaches to the investigation of vacuum (free) energy when the eigen- 
value problem is well-posed. In the first, one drops the zero-mode contribution. In usual 
quantum field theory this can be justified by the renormalization of the corresponding cou- 
pling constant (normally, the cosmological constant). In the second, one keeps the zero-mode 
contribution explicitly. Clearly, both these approaches are merely conventions and one must 
supply some physical considerations to argue in favor of one or the other. Moreover, the 
interpretation in terms of renormalization of coupling constants in quantum field theory at 
nonzero temperature becomes quite unclear, as all ultraviolet divergences are not influenced 
by temperature. 

In the present paper we adopted the first point of view and dropped the zero-mode 
contribution, which is consistent also with earlier calculations by Kutasov and Larsen ||, 
as well as Klemm, Petkou and Siopsis || |i~3| . After this paper appeared in hep-th, this 
approach was criticized by Dowker in Ref . |21[ , who argued in favor of the second convention, 
of keeping the zero-mode contribution. Dowker says that for the d = 2 conformal scalar, 
instead of Eq. ( |3.2| ) , we should have (there is an overall factor of 2 difference because he uses 
parity-degenerate states) 

U = -— + -Y'-fi£ , (Cl) 



12a 



n=0 



so there is an extra T term at low temperature which of course cancels the — T term at 
high temperature. [See Eq. (|3.3|) .1 Of course his remark about the Poisson sum formula 
being the basic mechanism does not establish the existence of the zero-mode contribution; 
rather, it just relates the two limits. Let us try to determine which convention is preferable 
by thermodynamical considerations. 

The first question is, where does the n = term come from? Of course there is a zero 
mode for a circle. But its contribution to the partition function seems ambiguous: 



InZ 



zero mode 



In (1 



-130 



) = - In /3 ■ 0, 



(C2) 
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so if one does not worry about the meaning of this, one gets 



d 1 

-Ezeromodc = TTTT ^ P = ~ni (^3) 

op p 



as Dowker states. This does not seem very well-defined. We would interpret Eq. ( |C2|) as an 
infinite constant, which should be dealt with by the process of renormalization. Furthermore, 
in Sec. Ill, we presented an independent derivation of the energy starting from the Green's 
function which directly gave the high-temperature limit ( |3.4| ) with the — T there. It seems 
clear there is a certain level of ambiguity here. 

This ambiguity is really only a mathematical one. Physically the situation is clear-cut: 
on a circle the eigenvalues are \n\/a, and thus zero for n = 0. This means that n = cannot 
contribute to the sum at all; this case is eliminated physically from the start. 

To completely resolve this ambiguity we note that such a term is in conflict with funda- 
mental thermodynamical principles. If for low temperature 

E = U-E = T= ^L/3F, (C4) 

the free energy must be, consistent with Eq. (|C2]), 

F = i (In 13 + constant) , (C5) 

P 

which then means that the entropy is 

S = [3(E - F) = 1 - In p - constant, (C6) 

which diverges to — oo as (3 — > oo. This grossly violates the Third Law of thermodynamics. 
In contrast, our temperature correction vanishes exponentially in that limit. [See Eq. (|6.8cj) .1 
Dowker bases his analysis on his earlier Ref . |22| . There indeed he obtains the free energy 



in the form ( |C5|) as a singular limit — i.e., the constant is infinite. He further mentions the 
entropy, saying that dropping the divergence is just a shift of the zero of entropy and hence 
is unobservable, and then stating that "it is not clear whether or not there is a problem with 
Nernst's theorem." Evidently there is such a problem, and appealing to the nonexistence of 
an ideal gas state seems quite beside the point. 

After this paper was completed, a preprint |2B| appeared which completely supports our 
conclusions. 
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